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Abstract 

In this review paper we discuss the effects of pairing correlations on inner crust 
matter in the density region where nuclear clusters are supposed to coexist with non- 
localised neutrons. The pairing correlations are treated in the framework of the finite 
temperature Hartree-Fock-Bogoliubov approach and using zero range nuclear forces. 
After a short introduction and presentation of the formalism we discuss how the pairing 
correlations affect the structure of the inner crust matter, i.e., the proton to neutron 
ratio and the size of Wigner-Seitz cells. Then we show how the pairing correlations 
influence, though the specific heat of neutrons, the thermalization of the crust in the 
case of a rapid cooling scenario. 
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1. Introduction 



The superfluid properties of the inner crust of neutron stars have been considered long ago 
in connection to the large relaxation times which follow the giant glitches. Thus, according 
to the present models, the glitches are supposed to be generated by the unpinning of the 
superfluid vortex lines from the nuclear clusters immersed in the inner crust of neutron 
stars H]|47l. Later on the superfluidity of the inner crust matter was also considered in 
relation to the cooling of isolated neutron stars 041 [T2l and, more recently, in the thermal 
after-burst relaxation of neutron stars from X-ray transients ll56l[T3ll2"9Tl . 

The superfluid properties of the inner crust are essentially determined by the non- 
localized neutrons. For baryonic densities smaller than about 1.4 x 10 -14 g cm -3 the 
non-localized neutrons are supposed to coexist with nuclei-type clusters IPTT1I431 . At higher 
densities, before the nuclear matter becomes uniform, the neutrons and the protons can form 
other configurations such as rods, plates, tubes and bubbles ll46l . 

A microscopic ab initio calculation of pairing in inner crust matter should take into 
account the polarization effects induced by the nuclear medium upon the bare nucleon- 
nucleon interaction. This is a very difficult task which is not yet completely solved even for 
the infinite neutron matter. Thus, compared to BCS calculations with bare nucleon-nucleon 
forces, most of variational or diagrammatic models predict for infinite matter a substantial 
reduction of the pairing correlations due to the in-medium polarisation effects (36). On the 
other hand, calculations based on Monte Carlo techniques predict for dilute neutron matter 
results closer to the BCS calculations (for a recent study see 11231 ). 

A consistent treatment of polarization effects on pairing is still missing for inner crust 
matter (for a recent exploratory study see [6]). Therefore at present the most advanced mi- 
croscopic model applied to inner crust matter remains the Hartree-Fock-Bogoliubov (HFB) 
approach. Pairing correlations have been also considered in the Quasiparticle Random 
Phase Approximation (QRPA)(see Section 2.1 below) in relation to the collective modes 
in inner crust matter OTTl . However, a systematic investigation of the effect of collective 
QRPA excitations on thermodynamic properties of inner crust matter is still missing. 

The scope of this chapter is to show how the HFB approach can be used to investigate 
the effects of pairing correlations on inner crust matter properties. Hence, in the first part 
of the chapter we will discuss the influence of pairing, treated in HFB approach at zero 
temperature, on the structure of inner crust matter. Then, using the HFB approach at finite 
temperature, we will show how the pairing correlations affect the specific heat and the 
thermalization of the inner crust matter in the case of a rapid cooling scenario. 

In the present study we will focus only to the region of the inner crust which is sup- 
posed to be formed by a bbc crystal lattice of nuclear clusters embedded in non-localized 
neutrons. The crystal lattice is divided in elementary cells which are treated in the Wigner- 
Seitz approximation. 
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2. Treatment of pairing in the inner crust of neutron stars 



2.1. Finite-temperature Hartree-Fock-Bogoliubov approach 

In this section we discuss the finite-temperature HFB approximation for a Wigner-Seitz cell 
which contains in its center a nuclear cluster surrounded by a neutron gas. The cell contains 
also relativistic electrons which are considered uniformly distributed. 

In principle, the HFB equations should be solved by respecting the bbc symmetry of the 
inner crust lattice. However, imposing the exact lattice symmetry in microscopic models 
is a very difficult task (for approximative solutions to this problem see Refs. lfl~8l l26l and 
the references therein). We therefore solve the HFB equations for a spherical WS cell, 
as commonly done in inner crust studies fi3l l3l. Since we are interested to describe the 
thermodynamic properties of the inner crust matter, we present here the HFB approach at 
finite temperature. 

The HFB equations for a spherical WS cells have the same form as for isolated atomic 
nuclei. Thus, for zero range pairing forces and spherical symmetry, the HFB equations at 
finite temperature are defined as |[24ll . 



h T , q {r) - \ 



Ar )9 (r) 



At, 9 M -h T , q (r) + X q 




(1) 



where E i>q is the quasiparticle energy, E i>q = J (e i)9 - A^) 2 + A? 9 , U i>q (r) and V itq (r) are 
the components of the HFB wave function and X q is the chemical potential ( q = n, p is the 
index for neutrons and protons) . The quantity hr,q[r) is the thermal averaged mean field 
hamiltonian and Ay^(r) is the thermal averaged pairing field. 

In a self-consistent HFB calculation based on a Skyrme-type force, as used in the 
present study, hr, q {r) and Aj. 5 (r) are expressed in terms of thermal averaged densities, 
i.e., particle density pT,q{r), kinetic energy density Tr,q{r), spin density Jr,q{r) and, re- 
spectively, pairing density KT, q {r). The thermal averaged densities mentioned above are 
given by 11531 : 
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where f,^ q = [1 + exp(Ei q /T)] 1 is the Fermi-Dirac distribution of quasiparticles, T is the 
temperature expressed in energy units, and gi q = 2ji tq + 1 is the degeneracy of the state i 
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with angular momentum jj. The summations in the equations above are over the spectrum 
of bound nucleons which form the nuclear cluster and of unbound neutrons which form the 
neutron gas. A constant density at the edge of the WS cell is obtained imposing Dirichlet- 
Von Neumann boundary conditions at the edge of the cell (43J, i.e., all wave functions of 
even parity vanish and the derivatives of odd-parity wave functions vanish. 

The nuclear mean field has the same expression in terms of densities as in finite nu- 
clei fl9l . However, for a WS cell the Coulomb mean field of protons has an additional 
contribution coming from the interaction of the protons with the electrons given by 

<Ur) = -e 2 Jd*r' Pe (r>) ] -L Tl . (6) 

Assuming that the electrons are uniformly distributed inside the cell, with the density p e = 

3Z/(4irRy VS ), one gets 

«SUO = - w* (nls - ^) = ^ ((^) 2 - 3) <7, 

It can be seen that inside the WS cell the contribution of the proton-electron interaction to 
the proton mean field is quadratic in the radial coordinate. 

The pairing field is calculated with a zero range force of the following form 

Vpair, g (ri, r-j) = V ffPair,g[PT,n(r),PTj)(r)](l - P<r)5(rij) , (8) 

where P a = (1 + a\ ■ 02) /2 is the spin exchange operator. For this interaction the pairing 
field is given by 

&T,q(r) = V 5Pair,g[PT,n(r),PT,p(r)] «r, 5 (r). (9) 

In the calculations presented here we use two different functionals for 
<7Pair,<j [pt,u (r), PT,p(^))- The first one, called below isoscalar (IS) pairing force, de- 
pends only on the total baryonic density, pt,b{t) = pr.n{f) + PT,p{r)- Its expression is 
given by 

3Pair>T,n(r),PT, P (r)] = 1 " V (^^f > ( 10 ) 

where po is the saturation density of the nuclear matter. This effective pairing interaction is 
extensively used in nuclear structure calculations and it was also employed for describing 
pairing correlations in the inner crust of neutron stars ll52l l53l l54l l40l . The parameters 
are chosen to reproduce in infinite neutron matter two pairing scenarii, i.e., corresponding 
to a maximum gap of about 3 MeV (strong pairing scenario, hereafter named ISS) and, 
respectively, to a maximum gap around 1 MeV (weak pairing scenario, called below ISW). 
These two pairing scenarii are simulated by two values of the pairing strength, i.e., Vo={- 
570,-430} MeV fm~ 3 . The other parameters are taken the same for the strong and the weak 
pairing, i.e., a=0.45, 7/=0.7 and po=0-16 fm -3 . The energy cut-off, necessary to cure the 
divergence associated to the zero range of the pairing force, is introduced through the factor 
g-Bi/ioo ac tj n g f or £j t > 20 MeV, where Ei are the HFB quasiparticle energies. 

The second pairing functional, referred below as isovector strong pairing (IVS), de- 
pends explicitly on neutron and proton densities and has the following form in the neutron 
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Figure 1. (color online) Neutron pairing gap for the interactions ISW (isoscalar weak), ISS 
(isoscalar strong) and IVS (isovector strong) in symmetric nuclear matter and in neutron 
matter. 



channel ll38l . 

0Pair,nK„(r),PT )P (r)H^ , (11) 

V Po J V Po J 

where I(r) = pT,n{r) — pT,p{r). This interaction is adjusted to reproduce the neutron 1 So 
pairing gap in neutron and symmetric nuclear matter provided by the BCS calculations with 
the bare nucleon-nucleon forces lfl4l . In addition, the pairing strength Vq and the cut-off 
energy are related to each other through the neutron-neutron scattering length according 
to the procedure described in Ref. 0. Therefore this interaction is expected to describe 
properly the pairing for all the nuclear densities of the inner crust matter, including the 
low density neutron gas. As shown in Refs. ll39l[T0ll . this pairing functional describes well 
the two-neutron separation energies and the odd-even mass differences in nuclei with open 
shells in neutrons. In the present calculations for this pairing functional we have used the 
parameters V = -703.86 MeV mi" 3 , r? s =0.7115, a s =0.3865, r? n =0.9727, a n =0.3906. The 
cut-off prescription is the same as for the isoscalar pairing force. 

The pairing gaps in symmetric matter and neutron matter predicted by the three pairing 
forces introduced above are represented in Fig.Q]for a wide range of sub-nuclear densities. 
It can be seen that the isovector IVS interaction gives a maximum gap closer to the strong 
isoscalar ISS force, and the ISW interaction predict a suppression of the pairing gap up to 
saturation density. 

To illustrate how the pairing correlations are spatially distributed in the Wigner-Seitz 
cells and how they are affected by the temperature, in Fig. 2 are shown the pairing fields 
for neutrons in the cells 2 and 5 (see Table 1). It can be noticed that the clusters have a non- 
trivial influence on the pairing of neutron gas. Thus, depending of the relative intensity of 
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Figure 2. Neutron pairing fields for the Wigner-Seitz cells 2 and 5 (see Table 1) labeled, 
respectively, as 1800 Sn and 950 Sn. The full and the long-dashed lines correspond to the 
ISS and ISW pairing interaction. Here the calculations have been done with the strenghts 
Vb = {-430, -330} MeV fm 3 and with an energy cut-off of 60 MeV. 



pairing in the cluster and the gas region, the presence of the cluster can suppress or enhance 
the pairing in the surface region of the cluster. 

2.2. Quasiparticle Random Phase Approximation (QRPA) 

Pairing correlations affect not only the ground state properties of inner crust matter but also 
its excitations modes. The non-collective excitations are commonly described by the quasi- 
particle energies obtained solving the HFB equations. To calculate the collective excitations 
one needs to take into account the residual interaction between the quasiparticles. In what 
follows we discuss briefly the collective modes of the inner crust matter in the framework 
of QRPA, which takes properly into account the pairing correlations OTTl . 

The QRPA can be obtained from the time-dependent HFB approach in the limit of linear 
response. In the linear response theory the fundamental quantity is the Green function 
which satisfies the Bethe-Salpeter equation 

G = (1-G V)- 1 G = G + G VG. (12) 

The unperturbed Green's function Go has the form: 

G «W^)=| 03) 

where E{ are the HFB quasiparticle energies and U%j are 3 by 2 matrices expressed in term of 
the two components of the HFB wave functions |[32l . The J] symbol in the equation above 
indicates that the summation is taken over the bound and unbound quasiparticle states. The 
latter corresponds here to the non-localised neutrons in the WS cell. 

V is the matrix of the residual interaction expressed in terms of the second derivatives 
of the HFB energy functional, namely: 
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Figure 3. Strength distributions for the quadrupole mode in WS cell 2 of Table 1. The 
full (dashed) line corresponds to the QRPA (HFB) strength. The results are for ISS pairing 
force with the same parameters as used in Figure 1. 



In the above equation {pi, p%, p%} = {p, k, k*}, where p and k are, respectively, the parti- 
cle (2) and pairing (5) densities; the notation a means that whenever a is 2 or 3 then a is 3 
or 2. 

The linear response of the system to external perturbation is commonly described by 
the strength function. Thus, when the external perturbation is induced by a particle-hole 
external field F the strength function writes: 

S(u) = --Im [ F*(r)G u (r,r';uj)F(r')dr dr' (15) 
vr J 

where G 11 is the (ph,ph) component of the QRPA Green's function. 

As an example in Fig. 3 it is shown the strength function for the quadrupole response 
calculated for the WS cell 2 of Table 1 below |f3TTl . The results correspond to the isoscalar 
pairing force with the strength Vo=-430 MeV fm~ 3 and an energy cut-off of 60 MeV. It 
can be seen that the unperturbed spectrum, distributed over a large energy region, becomes 
concentrated almost entirely in the peak located at about 3 MeV when the residual interac- 
tion between the quasiparticles is introduced. The peak collects more than 99% of the total 
quadrupole strength and it is extremely collective. An indication of the extreme collectivity 
of this low-energy mode can be also seen from its reduced transition probability, B(E2), 
which is equal to 25 x 10 3 Weisskopf units. This value of B(E2) is two orders of magnitude 
higher than in standard nuclei. This underlines the fact that in this WS cell the collective 
dynamics of the neutron gas dominates over the cluster contribution. In Ref. I 3H it is shown 
that similar collective modes appears for the monopole and the dipole excitations. A very 
collective low-energy quadrupole mode it was also found in all the Wigner-Seitz cells with 
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Table 1 . The structure of the Wigner-Seitz cells obtained from a density matrix expansion 
(DME) l43l pb is the baryon density, kp = (S^ne) 1 / 3 the Fermi momentum calculated 
as in Ref. [3] where ns is the number of baryons per fm 3 , N and Z are the numbers of 
neutrons and protons while Rws is the radius of the WS cells. Compared to Ref. [43 ] here 
it is not shown the cell with the highest density located at the interface with the pasta phase. 

Z=50 [27]. However, a systematic investigation of the influence of these collective modes 
on the thermodynamic properties of inner crust matter is still missing. 

3. The effect of pairing on inner crust structure 

The first microscopic calculation of the inner-crust structure was performed by Negele and 
Vautherin in 1973 [43]. In this work the crystal lattice is divided in spherical cells which 
are treated in the Wigner-Seitz (WS) approximation. The nuclear matter from each cell is 
described in the framework of Hartree-Fock (HF) and the pairing is neglected. The proper- 
ties of the WS cells found in Ref. |43l . determined for a limited set of densities, are shown 
in Table I. The most remarkable result of this calculation is that the majority of the cells 
have semi-magic and magic proton numbers, i.e., Z=40,50. This indicates that in these 
calculations there are strong proton shell effects, as in isolated atomic nuclei. 

The effect of pairing correlations on the structure of Wigner-Seitz cells was first inves- 
tigated in Refs. El|4l|5l within the Hartree-Fock BCS (HFBCS) approach. In this section 
we shall discuss the results of a recent calculations based on Hartree-Fock-Bogoliubov ap- 
proach 11281 . This approach offers better grounds than HF+BCS approximation for treating 
pairing correlations in non-uniform nuclear matter with both bound and unbound neutrons. 

As in Ref. [43 ], the lattice structure of the inner crust is described as a set of independent 
cells of spherical symmetry treated in the WS approximation. For baryonic densities below 
p m 1.4 x 10 14 g/cm 3 , each cell has in its center a nuclear cluster (bound protons and 
neutrons) surrounded by low-density and delocalized neutrons and immersed in a uniform 
gas of ultra-relativistic electrons which assure the charge neutrality. At a given baryonic 
density the structure of the cell, i.e., the N/Z ratio and the cell radius is determined from 
the minimization over N and Z of the total energy under the condition of beta equilibrium. 
The energy of the cell, relevant for determining the cell structure, has contributions from 
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the nuclear and the Coulomb interactions. Its expression is written in the following form 



E = Em + En + T e + El- 



(16) 



The first term is the mass difference Em = Z(m p + m e ) + (N — A)m n where N and Z are 
the number of neutrons and protons in the cell and A=N+Z. En is the binding energy of 
the nucleons, which includes the contribution of proton-proton Coulomb interaction inside 
the nuclear cluster. T e is the kinetic energy of the electrons while El is the lattice energy 
which takes into account the electron-electron and electron-proton interactions. The contri- 
bution to the total energy coming from the interaction between the WS cells Pfll it is not 
considered since it is very small compared to the other terms of Eq.(l). 

We shall now discuss the effect of pairing correlations on the structure of the WS cells. 
To study the influence of pairing correlations we have performed HFB calculations with the 
three pairing interactions introduced in Section 2.1. 
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Figure 4. (color online) The HFB energies per particle as function of proton number for the 
pairing forces ISW (dotted line), ISS (dashed line) and IVS (dashed-dotted line). The solid 
lines represent the HF results. In the left pannel are shown the results obtained including 
the finite size corrections. 



The structure of the WS cells obtained in the HFB approach is given in Table II while 
in Fig. @] it is shown the dependence of the binding energies, at beta equilibrium, on protons 
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number. The contribution of pairing energy, nuclear energy and electron kinetic energy to 
the total energy are given in Fig. 5 for the cells 2 and 6. From this figure we can notice 
that the nuclear binding energy is almost compensated by the kinetic energy of the electrons 
which explains the weak dependence of the total energy on Z seen in Fig. 4. In Fig.4 we 
observe also that pairing is smoothing significantly the variation of the HF energy with Z. 
For this reason the HFB absolute minima are very little pronounced compared to the other 
neighboring energy values. 
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Table 2. The structure of Wigner-Seitz cells obtained in the HF and HFB approximation. 
The results corresponds to the isoscalar weak (ISW), isoscalar strong (ISS) and isovector- 
isoscalar (IVS) pairing forces. In the last 4 columns are shown the proton numbers obtained 
with the finite size corrections. In the table are shown only the structures of the cells which 
could be well-defined by the present calculations. 



From Fig. [4] we observe that in the cells 1-2 the binding energy does not converge to a 
minimum at low values of Z. For the cells 3-4, although absolute minima can be found for 
HF or/and HFB calculations, these minima are very close to the value of binding energy at 
the lowest values of Z we could explore. Therefore the structure of these cells is ambiguous. 
The situation is different in the cells 5-10 where the binding energies converge to absolute 
minima which are well below the energies of the configurations with the lowest Z values. 

From Table II it can be observed that the structure of some cells becomes very differ- 
ent when the pairing is included. Moreover, these differences depend significantly on the 
intensity of pairing (see the results for ISW and ISS). However, as seen in Fig. 4, for the 
majority of cells the absolute minima are very little pronounced relative to the other energy 
values, especially for the HFB calculations. Therefore one cannot draw clear quantitative 
conclusions on how much pairing is changing the proton fraction in the cells. 

When the radius of the cell becomes too small the boundary conditions imposed at the 
cell border through the WS approximation generate an artificial large distance between the 
energy levels of the nonlocalized neutrons. Consequently, the binding energy of the neutron 
gas is significantly underestimated. An estimation of how large could be the errors in the 
binding energy induced by the WS approximation can be obtained from the quantity 

f(Pn,R\Vs) = Bi n fXPn) ~ B WS -inf.(Pn, R\Vs) > ( l7 ) 

where the first term is the binding energy per neutron for infinite neutron matter of density 
p n and the second term is the binding energy of neutron matter with the same density 
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Figure 5. (color online) The different contributions to the total energy in the cells 2 and 6 
for the HF (left pannel) and HFB (right pannel) calculations. Are shown: the total energy 
(solid line), the nuclear energy (dashed line), the kinetic energy of the electrons (dashed- 
dotted line), and the pairing energy for the ISS pairing interaction (dotted line). The pairing 
energies are shifted up as indicated in the figure. 



calculated inside the cell of radius Rws an d employing the same boundary conditions as in 
HF or HFB calculations. In Ref. 11371 it was proposed for the finte size energy correction, 
Eq. 26, the following parametrisation 

f(Pn a ,Rws) = 89.05(p ng /p f U25 Rw 2 s , (18) 

where p ng is the average density of neutrons in the gas region extracted from a calculation 
in which the cell contains both the nuclear cluster and the nonlocalized neutrons while po 
is the nuclear matter saturation density. 

How the energy corrections described by Eq. (18) influences the HF (HFB) results 
can be seen in Fig. 4 (right pannel) and in Table II (last four columns). As expected, the 
influence of the corrections is more important for the cells 1-5, in which the neutron gas has 
a higher density, and for those configurations corresponding to small cell radii 



4. Effect of pairing on the thermodynamic properties in the in- 
ner crust matter 

4.1. Specific heat of the baryonic matter 

The specific heat of the inner crust has contributions from the electrons, the lattice (or Ion) 
and the baryonic matter (essentially the de-localized neutrons). The contribution of these 
3 components is shown in Fig. [6] where different scenarios for the neutrons are considered. 
If neutrons are non-superfluid they dominate the contribution of the other component in the 
inner crust. The effect of the neutron superfluidity is the reduce strongly the contribution of 
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Figure 6. Specific heat in the crust of the neutron star for the different components of the 
star matter and for a temperature of 10 9 K. 



the neutrons. It is therefore expected a very large effect induced by the neutron superfluidity 
that we now discuss. The specific heat of the neutrons is calculated from the quasiparticle 
energies obtained solving the HFB equations at finite-temperature in the Wigner-Seitz ap- 
proximation, as discussed in Section 2.1. For the WS cells we use the structure determined 
inRef.Ea (see Table 1). 

The specific heat of the neutrons is calculated by 

where V is the volume of the Wigner-Seitz cell, T is the temperature and S is the entropy. 
The latter is given by 

S = -k B J2 9iMi,<i m ki + C 1 - hi) ln (! " hi))- ( 20 ) 

i,q 

In Fig. 7 are shown the results obtained for the specific heat in the case of a strong 
isoscalar pairing interaction. To illustrate the particular behavior of the specific heat in 
non-uniform matter and the validity of various approximations, in what follows we shall 
discuss in more detail the results for the cell number 6, which contains N=460 neutrons 
and Z=40 protons (see Table 1). In this cell the HFB calculations predict 378 unbound 
neutrons. The specific heat given by the HFB spectrum, in which the contribution of the 
cluster is included, is shown in Fig. [8] by full line. In the same figure are shown also 
the specific heats corresponding to two approximations employed in some studies 1341 l48l 
l45l . In these approximations the non-uniform distribution of neutrons is replaced with 
a uniform gas formed either by the total number of neutrons in the cell (N=460, dashed 
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Figure 7. Specific heat of neutrons in various WS cells for isoscalar strong ISS pairing. The 
specific heat is given in units of Boltzmann constant ks- 
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Figure 8. Neutron-specific heat in the the WS cell 6 for isoscalar strong pairing. The results 
corresponds to various approximations discussed in the paper. The specific heat is given in 
units of Boltzman constant ks- 



line), or by taking only the number of the unbound neutrons (N=378, dashed-dotted line). 
How these approximations work is seen in Fig. [8] To make the comparison meaningful, 
the calculations for the uniform neutron gas are done solving the HFB equations with the 
same boundary conditions as for the non-uniform system, i.e., neutrons+cluster. As seen 
in Fig. [8j the transition from the superfiuid phase to the normal phase is taking place at a 
lower temperature in the case of uniform neutron gas, especially when are considered only 
the unbound neutrons. The critical temperature is therefore lower in uniform matter (for the 
two prescriptions often used) than in non-uniform matter. 

4.2. Pairing and thermalization time of the inner crust 

Several studies have shown that the thermalization time depends significantly on the super- 
fluid properties of the inner crust baryonic matter ll34l l25l l48l l40l l54ll . This dependence 
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Figure 9. Time evolution of the apparent surface temperature for the initial temperatures 
Ti = 500 keV. NC indicates the results of the calculations obtained by neglecting the effect 
of the clusters. 



is induced through the specific heat of unbound neutrons, strongly affected by the pairing 
energy gap. As an example we refer to Fig. 6 where it is seen that neutron superfluid- 
ity strongly suppress the neutron specific heat which, in the absence of the pairing would 
dominate over the specific heat of electrons and the lattice. 

To estimate the effect of neutron superfluidity to the thermalization time of inner crust 
matter we use a rapid cooling scenario in which the core arrives quickly to a much smaller 
temperature than the crust due to the direct URCA colling process. The thermalization 
time is defined as the time for the core nd the crust temperatures to equilibrate. The heat 
diffusion is described by the relativistic heat equation ll57l : 



d_ 

Or 



T(r) dr 



r 2 r(r)e* (c v ^ + e^X (21) 



where t is the time, K is the thermal conductivity, Cy is the specific heat and Q v is the 
neutrino emissivity. The effect of the gravity is given through the gravitational potential 
(f), which enters in the definition of the redshifted temperature T = Te^, and the quantity 
T(r) = (1 — 2Gm(r)/rc 2 ) 1 ^ 2 , where G is the gravitational constant and m(r) is the 
gravitational mass included in a sphere of radius r. The latter is obtained from the Tolman- 
Oppenheimer-Volkoff (TOV) equations [29] based on an equation of state obtained from 
SLy4 Skyrme interaction GUI . More details on the cooling model are given in Ref. Il22ll . 

The time evolution of the apparent surface temperature T~ is displayed in Fig. [9] for 
the initial temperatures of the crust T, = 500 keV ll22l . The effective surface temperature 
shown in Fig. [9] is obtained from the temperature at the bottom of the crust, Tb = T(p^), 
where p\ y = 10 10 g.cm -3 , using the relationship given in Ref. [50] for a non-accreted 
envelope. The results shown in Fig. [9] corresponds to a neutron star of mass 1.6 M Q in 
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Figure 10. Cooling times i w versus scaling parameter a for three pairing scenarii as dis- 
cussed in the text. The results correspond to neutron stars with masses between 1.4 and 2.0 
M & . The fitting curves are given for the case of a linear scale (dashed lines, right side) and 
for a fractional power of the scaling parameter a (solid lines, left side). 



which the inner crust extends from R c = 10.72 km, which is the radius at the core-crust 
interface, to 11.19 km. As can be noticed, the pairing enhances significantly the cooling at 
the surface of the star. In Fig. [9] are also shown the apparent surface temperatures obtained 
neglecting the effect of the clusters, i.e., supposing that the neutron specific heat is given 
solely by that of the neutron gas. In this case the neutron specific heat is calculated from 
the quasiparticle spectrum of BCS equations solved for infinite neutron matter at a density 
corresponding to that of the external neutrons in the WS cell [35]. In the case of weakly 
pairing scenario (ISW), the apparent surface temperature is dropping faster for superfluid 
non uniform matter than for superfluid uniform matter. For the strong pairing scenario, 
since the pairing correlations suppress the role of the neutron in almost the entire inner 
crust, the effect of the clusters is less important [22]. 

A simple random walk model for the cooling is sometimes used to estimate the thermal- 
ization time l[T2ll48ll40l . In this model the diffusion of the cooling wave towards the core 
is calculated without taking into account the dynamical change of the temperature through 
the whole crust and neglecting the neutrino emissivity in Eq. (121) . The crust is divided into 
shells of thickness Ri for which the thermal diffusively Di = Ki/Cy,i is considered as 
constant. From a dimensional analysis of the heat equation (|2TI) . the relaxation time of each 
of these shells is defined as, 

n = r( J R) 3 J R l 2 /( 7 A) ! (22) 

where 7 is a geometrical factor and R is the radius of the neutron star. The total relaxation 
through the crust is given by, 
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The thermalization time defined as d23l satisfy the condition that in a uniform density case, 
Tth is independent of the shell discretization. Using this model, the effect of the neutron 
superfluidity has also been estimated to be large H8ll40l . and is qualitatively similar to the 
result deduced from solving of the heat diffusion equation (f2TT >. 

In some previous studies IT341 l25l it was found that the thermalization time, defined as 
the time needed for the crust to arrive to the same temperature with the core, satisfies the 
scaling relation = t\oP where 



depends solely on the global properties of the neutron star, i.e., the crust thickness Ai? cmst , 
the star radius R and the mass of the star M. The scaling relation can be inferred from the 
random walk model presented above. In Fig. \\0\ are shown as a function of the scaling 
parameter a the cooling times we have obtained from the solution of the heat diffusion 
equation (1211 ) for various neutron stars masses lying between 1.4 and 2.0 M G with a step 
of O.IMq. In this figure are displayed two sets of fits for t w versus a, i.e., a linear fit 
with (3 = 1 (dashed line) and a fit with a fractional value for (3 (solid line). It can be 
noticed that the best fit is obtained with a fractional value for j3, which is equal to 0.86 
for the normal neutrons, to 0.85 for weakly paired neutrons and to 0.89 for strongly paired 
neutrons. Considering the simpler linear fit, i.e., t w rj od\, as done in 04ll25l . we get for 
the normalized time t\ the values t\ = {68.9,39.3,22.3} corresponding, respectively, to 
the normal neutrons, neutrons with weak pairing and neutrons with strong pairing. For a 
1.5 M & neutron star with a = 1.15, we obtain t\ = {66.4, 37.4, 26.6}. These values for t\ 
are larger than that of Ref. l25l (Table 2) by a factor 2.3 in the non-superfluid case, and 3.4 
(3.0) for the weak (strong) pairing scenario. These differences could be explained by the 
effects of the nuclear clusters on the neutron specific heat, disregarded in Ref. 11251 . and by 
different neutrinos processes and thermal conductivities in the core matter used in the two 
calculations. 

5. Summary and Conclusions 

In the first part of this chapter we have discussed the influence of pairing correlations on 
the structure of inner crust of neutron stars. The study was done for the region of the 
inner crust which is supposed to be formed by a lattice of spherical clusters embedded in 
a gas of neutrons. The lattice was treated as a set of independent cells described in the 
Wigner-Seitz approximation. To determine the structure of a cell we have used the nuclear 
binding energy given by the Skyrme-HFB approach. For the cells with high density and 
small radii the binding energies do not converge to a minimum when the proton number 
has small values. We believe that it is related to the discretization of the continuum. For a 
small radius of the cell the average distance between the energy levels of the non-localised 
neutrons becomes artificially large which cause an underestimation of the binding energy. 
To correct this drawback we have used an empirical expression based on the comparison 
between the binding energy of neutrons calculated in infinite matter and in a finite-size 
spherical cell [37]. We found that the finite size corrections to the binding energies are 




(24) 
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significant for the high density cells with small proton numbers. This conclusion indicates 
the need of a more accurate evaluation of the errors induced by the finite size of the WS 
cells on nuclear binding energy, which requires to go beyond the empirical expression used 
here. 

In the second part we have discussed how thermalization of neutron stars crust depends 
on pairing properties and on cluster structure of the inner crust matter. The thermal evolu- 
tion was obtained by solving the relativistic heat equation with initial conditions specific to 
a rapid cooling process. The specific heat of neutrons was calculated from the HFB spec- 
trum. The results show that the crust thermalization is strongly influenced by the intensity 
of pairing correlation. It is also shown that the cluster structure of the inner crust affects 
significantly the time evolution of the surface temperature. 

The thermodynamic properties of the inner crust discussed in this chapter are based 
on the excitation spectrum of HFB equations. As it was shown in Section 2.2, the QRPA 
approach predicts in some cells collective excitations located at low energies, comparable 
with the pairing energy gap. It is thus expected a significant contribution of the collective 
modes to the specific heat and to the thermalization process of inner crust matter. This is an 
open issue which deserves further investigations. 
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